Device and system for detecting abnormality

ABSTRACT

An abnormality detection device  30  includes small motion sensors  25   a˜   25   c  that detect small motions of a person in a house, a data collecting unit  32  that collects and stores sensor signals from the small motion sensors  25   a˜   25   c  as sensor patterns, a Markov chain operating unit  33  that transforms the sensor patterns into a cluster sequence by vector-quantizing input patterns which are obtained by averaging and normalizing the sensor patterns and calculates a transition number matrix and a duration time distribution of a Markov chain and so on using a Markov chain model, a comparing unit  34  that calculates characteristic amount (Euclid distance and average log likelihood in appearance frequency of a Markov chain and average log likelihood to the duration time distribution of a Markov chain) of a sample activity as against a daily activity based on the obtained transition number matrix and the duration time distribution and so on, and others.

BACKGROUND OF THE INVENTION

[0001] (1) Field of the Invention

[0002] The present invention relates to an abnormality detection device, particularly to a device and others suitable for detecting abnormality such as an unusual human activity which has occurred in an elderly person who lives alone in a house.

[0003] (2) Description of the Related Art

[0004] Various attempts have been made for detecting abnormality which has occurred in an event such as machine operation and human behavior. For example, it is an important technology, particularly in rapidly aging Japan, to monitor activities of an elderly person who lives alone and detect occurrence of his unusual activity. One of such conventional technologies is “Abnormality Report System” disclosed in the Japanese Laid-Open Patent Application No. 2001-67576.

[0005] This conventional system includes, for monitoring the daily life of an elderly person or a sick person who lives alone, (a) a sensor unit that is placed in the restroom in the house of the person subject to monitoring and outputs a predetermined signal when detecting the use of the restroom, (b) a first communication means (a wireless terminal device) that is placed in the restroom and outputs a predetermined signal when receiving the signal from the sensor unit, and (c) a second communication means (a main device) that includes a monitoring timer which starts clocking when receiving the signal from the first communication means and reports the occurrence of abnormality to the monitoring center when the signal from the first communication means is interrupted for a predetermined time period or more. When the elderly person faints due to a disease or the like and cannot move at all, he does not use the restroom for the predetermined time period or more, and thereby, it is detected that something abnormal has occurred in him.

[0006] However, the above-cited conventional system is based on his behavior characteristic that he always goes to the restroom within the predetermined time period. Therefore, it requires some steps to eliminate an exceptional case when the elderly person does not go to the restroom for the predetermined time period or more even in nothing abnormal (normal situation), such as the cases when he goes out and he does his needs with a chamber pot.

[0007] Furthermore, since this conventional system only detects unusual human activities in a building with a restroom, there is a problem that the uses and available opportunities thereof are extremely limited. In other words, the conventional system can detect neither unusual human activities in a place without a restroom or in a time period when the restroom is not used, nor abnormal motions of animals such as a pet or machines which do not use a restroom.

SUMMARY OF THE INVENTION

[0008] In order to solve the aforesaid problem, the object of the present invention is to provide an abnormality detection device and others with a variety of applications that can detect occurrence of abnormality of objects and events under observation, independently of the types and numbers thereof and the space (such as a place and a time zone) where they are observed.

[0009] In order to achieve the object, the abnormality detection device according to the present invention is a device that detects occurrence of abnormality in an event under observation, comprising: an input unit operable to acquire a sequence of an input pattern that is data which depends upon the event; a transition analyzing unit operable to analyze a characteristic of a transition in the acquired sequence of the input pattern; a comparing unit operable to compare the analyzed characteristic of the transition with a predetermined reference value, and judge that abnormality has occurred in the event when the characteristic and the reference value are not approximate to each other within a predetermined range; and an output unit operable to output occurrence of abnormality when the comparing unit judges that abnormality has occurred.

[0010] In other words, the occurrence of abnormality is detected by acquiring a sequence (time series) of input patterns that are the data depending upon changes of an event such as a motion of a person or a thing to be observed, focusing on the transition of the patterns in the sequence, sampling the characteristic amount in the transition and comparing it with that in the normal case. Here, the “event” means an event that can be represented as data which can be processed by a computer based on a signal from an equipment sensor, a report (data entry) by a person, and others, and is typically a “human activity” in a house.

[0011] When abnormality is detected based on the human activities in the house, the following approach is taken in the present invention.

[0012] Description of a personal daily life and detection of unusual events in his life are made using sensor information obtained in a house (an intelligent house) where small motion sensors are placed in a plurality of predetermined places. The objective here is that a computer understands the daily life which is individually customized because human behavior varies greatly from person to person.

[0013] On the other hand, there is symbol processing where discretized environments are used for prediction. It is necessary for us to predict events in an outside world for our lives, and such an intellectual function may be used for modeling and information processing in the equipment which matches with human beings and supports them.

[0014] The intelligent house equipped with sensors that detect human presence and activities is considered to be something equipped with human outward sensory organs inwardly, and shows effectiveness as a place where the intellectual information processing for supporting human activities is applied. Therefore, based on the data actually measured in the intelligent house, unusual behavior is detected using the likelihood of the activity sequence as an information processing model similar to that of human beings.

[0015] More specifically, discretization of continuous sensory data that is a basis of symbol processing in human brains is performed, an environmental model is constructed based on the transition between the discrete states, and then the likelihood of the actual activity is predicted based on that model. Accordingly, an abnormal activity (an activity with low likelihood) can be found based on the comparison with the learned personal data, and therefore an abnormal activity may be detected without disturbing his daily life. It is preferable to use vector quantization for discretization and a Markov process model for the environmental model.

[0016] Furthermore, on the assumption that a human life is made up of a plurality of daily activities and they are triggered according to the situation, the sensor sequence is discretized into a reproducible time interval and a local daily activity template in that time interval is generated. Since an image processing method is applied to automatic sampling of time intervals, highly correlative time intervals are sampled by template matching using a time window. Thereby, daily life representation with a hierarchy in the time domain can be constructed.

[0017] As a verification of the daily activity template, unusual activities were detected. The unusual activities were detected by detecting the differences between the daily activity template and the actual activities. More specifically, methods such as (i) comparison between the probability distribution based on the daily activity template and the likelihood of the actual activities and (ii) measurement of the difference between the daily activity template calculated in the local time interval and the global daily activity template using the probability distribution distance are used so as to evaluate them.

[0018] The present invention is not only realized by dedicated hardware such as the above-mentioned abnormality detection device, but also realized as an abnormality detection method including steps of the characteristic constituent elements, or as a program that causes a general-purpose computer to execute these steps, or as an abnormality detection system including the abnormality detection device and the receiver of the abnormality report.

FURTHER INFORMATION ABOUT TECHNICAL BACKGROUND TO THIS APPLICATION

[0019] The following applications are incorporated herein by reference:

[0020] Japanese Patent Application Ser. No. 2001-392921 filed Dec. 25, 2001;

[0021] Japanese Patent Application Ser. No. 2002-111292 filed Apr. 12, 2002.

BRIEF DESCRIPTION OF THE DRAWINGS

[0022] These and other objects, advantages and features of the invention will become apparent from the following description thereof taken in conjunction with the accompanying drawings that illustrate a specific embodiment of the invention. In the Drawings:

[0023]FIG. 1 is a block diagram showing an overall structure of the abnormality detection system according to the embodiment of the present invention.

[0024]FIG. 2 is a plan view of an intelligent house showing locations of small motion detection sensors.

[0025]FIG. 3 is a functional block diagram showing a structure of the abnormality detection device.

[0026]FIG. 4 shows values of sensor signals from the small motion detection sensors (sequences of sensor patterns).

[0027]FIG. 5 is a flowchart showing operation of a Markov chain operating unit.

[0028]FIG. 6 is a flowchart showing a detailed calculation procedure for averaging sensor patterns.

[0029]FIG. 7 is a flowchart showing a detailed calculation procedure for normalizing the averaged sensor pattern.

[0030]FIG. 8 is a flowchart showing a detailed calculation procedure for vector quantizing input patterns.

[0031]FIG. 9 is a plan view of the intelligent house showing locations of clusters.

[0032]FIG. 10 is a flowchart showing a calculation procedure of a cluster sequence.

[0033]FIG. 11 is a flowchart showing a calculation procedure of a transition number matrix in appearance frequency of Markov chain.

[0034]FIG. 12 is a flowchart showing a first half of a calculation procedure of duration time distribution of Markov chain.

[0035]FIG. 13 is a flowchart showing the second half of the calculation procedure of the duration time distribution of the Markov chain.

[0036]FIG. 14 is a data flow diagram showing a calculation procedure of distance in the appearance frequency of Markov chain performed by the distance calculating unit of the comparing unit.

[0037]FIG. 15 is a flowchart showing a detailed calculation procedure for normalization in the above distance calculation.

[0038]FIG. 16 is a flowchart showing a detailed calculation procedure for calculating Euclid distance in the above distance calculation.

[0039]FIG. 17 is a data flow diagram showing a calculation procedure of likelihood in the appearance frequency of Markov chain performed by the likelihood calculating unit of the comparing unit.

[0040]FIG. 18 is a flowchart showing a detailed calculation procedure in the above likelihood calculation.

[0041]FIG. 19 is a data flow diagram showing a calculation procedure of likelihood for the duration time distribution of Markov chain performed by the likelihood calculating unit of the comparing unit.

[0042]FIG. 20 is a flowchart showing a detailed calculation procedure in the above likelihood calculation.

[0043]FIG. 21 is a plan view of the intelligent house showing the transition number matrix of a daily activity template in Experiment 1.

[0044]FIG. 22 is a plan view of the intelligent house showing the transition number matrix of an unusual activity in Experiment 1.

[0045]FIG. 23 is a diagram showing the distance of the appearance frequency of Markov chain for each experiment number (sample) in Experiment 1.

[0046]FIG. 24 is a diagram showing the distance of the appearance frequency of Markov chain for each experiment number (sample) in Experiment 2.

[0047]FIG. 25 is a diagram showing the likelihood for the duration time distribution of Markov chain in Experiment 3.

[0048]FIG. 26 is a plan view of the apartment style intelligent house showing the transition number matrix of the daily activity template in Experiment 4.

[0049]FIG. 27 is a diagram showing the likelihood of the appearance frequency of Markov chain on every experiment day and time in Experiment 4.

[0050]FIG. 28 is a diagram showing activities of the subject person based on his self record in Experiment 5.

[0051]FIG. 29 is a diagram showing the likelihood of the appearance frequency of Markov chain on every experiment day and time in Experiment 5.

[0052]FIG. 30 is a functional block diagram showing a detailed structure of the Markov chain operating unit of the abnormality detection system according to a modification of the present invention.

[0053]FIG. 31 is a flowchart showing operation of a local template generating unit of the abnormality detection system shown in FIG. 30.

[0054] FIGS. 32A˜D are diagrams showing small motion sensor sequences for explanation of a local template.

[0055] FIGS. 33A˜J are diagrams showing calculation results of correlation of all the times for 4 types of the sensor sequences as shown in FIGS. 32A˜D.

[0056] FIGS. 34A˜J are diagrams showing examples of binarization and labeling of the correlation shown in FIGS. 33A˜J and sampling of each time period thereof.

[0057] FIGS. 35A˜H are diagrams showing examples of sampling of similar sequence intervals from 4 types of the sequences for morning scenarios.

[0058] FIGS. 36A˜D are diagrams showing examples of sampling of similar sequence intervals from 2 types of the sequences for evening scenarios.

[0059]FIG. 37 is a flowchart showing operation of a local template predicting unit of the abnormality detection device shown in FIG. 30.

[0060] FIGS. 38A˜H are diagrams showing reliability and predicted templates.

[0061] FIGS. 39A˜H are diagrams showing input cluster sequences and predicted templates.

[0062] FIGS. 40A˜F are diagrams showing sequences for sampling 6 types of the local templates, FIGS. 40F˜I are diagrams showing local templates before selecting the local templates corresponding to FIGS. 40A˜F, and FIGS. 40M˜R are diagrams showing local templates after selecting the local templates corresponding to FIGS. 40A˜F.

[0063]FIG. 41 is a diagram showing a comparison of the average log likelihood between the case where a global template is used and a local template is used in Experiment 1.

[0064]FIG. 42 is a diagram showing a comparison of the average log likelihood between the case where the global template is used and the local template is used in Experiment 2.

[0065]FIG. 43 is a diagram showing a comparison of the average log likelihood between the case where the global template is used and the local template is used in Experiment 1/2 .

[0066]FIG. 44 is a diagram showing the configuration of the abnormality detection system according to a modification of the present embodiment that detects the occurrence of abnormality in human activities using various sensor signals.

[0067]FIG. 45 is a diagram showing the configuration of the abnormality detection system according to another modification of the present embodiment that performs decentralized processing in which the intelligent house collects data and the monitoring center or the like analyzes and judges the occurrence of abnormality.

[0068]FIG. 46 is a diagram showing the configuration of the abnormality detection system according to still another modification of the present embodiment including the abnormality report device that collects the sensor signals via the home network.

DESCRIPTION OF THE PREFERRED EMBODIMENT(S)

[0069] The embodiment of the present invention will be explained in detail with reference to the figures.

[0070]FIG. 1 is a block diagram showing an overall structure of the abnormality detection system according to the embodiment of the present invention. The abnormality detection system 10 is a system that detects occurrence of abnormality, that is, unusual human activities, in a house and report it to a specific contact station. The abnormality detection system 10 is characterized by detecting unusual activities without limiting the places in the house or the time of day, and includes an intelligent house 20, a monitoring center 50 and a cell phone 60 which are connected to each other via a communication network 40.

[0071] The intelligent house 20 is a two-story house shown in a plan view of FIG. 2, and includes small motion detection sensors 25 a˜25 c that are placed at predetermined locations (20 locations indicated by black circles) in the house and an abnormality detection device 30 that detects unusual human activities in the house by monitoring sensor signals outputted from the small motion detection sensors 25 a˜25 c and reports its detection results to the monitoring center 50 and a predetermined contact station (the cell phone 60 in this case) via the communication network 40. The small motion detection sensors 25 a˜25 c are infrared sensors that detect infrared rays emanating from the human body, for instance, and when there exists a person in the detection space (in the restroom, for instance) determined by the orientation of each sensor and they detect the motion of the person, they output it to the abnormality detection device 30.

[0072] The monitoring center 50 is a nursing-care center or the like that keeps centralized monitoring on occurrence of unusual activities in a plurality of the target intelligent houses 20, and includes a device that outputs to an alarm or a display device its receipt of the report from the abnormality detection device 30 via the communication network 40. The cell phone 60 is one of the contact stations registered in the abnormality detection device 30, and receives the report sent from the abnormality detection device 30 (such as an e-mail indicating that unusual activity has occurred in the intelligent house 20).

[0073]FIG. 3 is a functional block diagram showing the structure of the abnormality detection device 30 shown in FIG. 1. The abnormality detection device 30 is a controller which is connected to the small motion detection sensors 25 a˜25 c in the house and the communication network 40, and includes an input unit 31, a data collecting unit 32, a Markov chain operating unit 33, a comparing unit 34 and an output unit 35.

[0074] The input unit 31 is an input device such as an operation panel, and is used by the operator for giving the data to the abnormality detection device 30 instead of the small motion detection sensors 25 a˜25 c or setting various parameters for the abnormality detection device 30.

[0075] The data collecting unit 32 is comprised of a data logger or a hard disk that receives and records the sensor signals outputted from the small motion detection sensors 25 a˜25 c, and includes a sample activity data storage unit 32 a that stores the sensor signals outputted from the small motion detection sensors 25 a˜25 c and the input data outputted from the input unit 31 as the data of sample activities according to predetermined conditions, and a daily activity data storage unit 32 b that stores them as the data of daily activities.

[0076] Here, the “sample activity” means an activity that is to be a target of monitoring occurrence of abnormality, and the “daily activity” means an activity under the normal condition (the activity that is the measure (reference or standard) for judging whether unusual activity has occurred or not).

[0077] The Markov chain operating unit 33 is a processing unit which is realized by a CPU, a memory or the like that executes data processing based on the control program, and calculates the transition number matrix and the duration time distribution of Markov chain in the daily activities as characteristic amount by performing operations (which will be described later) on the sensor patterns stored in the daily activity data storage unit 32 b prior to monitoring the occurrence of unusual activity (in the mode of creating the reference data). In the mode of monitoring the occurrence of unusual activity (monitoring mode), the Markov chain operating unit 33 performs the same operation on the sequence (time series) of the sensor patterns stored in the sample activity data storage unit 32 a, or directly on the sensor patterns inputted by the small motion detection sensors 25 a˜25 c, so as to calculate the sequence of clusters, the transition number matrix of Markov chain and so on as the characteristic amount in the sample activities. The obtained results are respectively outputted to the comparing unit 34.

[0078] “Markov chain” means the sequence of the events when the nth event is determined in relation to the previous events in the sequence of the events, and in the present embodiment, it corresponds to the sequence of sensor pattern indicating the human activities (such as places where he moves) in the house (including the input pattern and the sequence of clusters which are obtained by performing a predetermined operation on that sequence). The “transition number matrix” of Markov chain is a matrix indicating the number of transitions from the past events to that event. The “duration time distribution” of Markov chain is a histogram of the duration time of each element (state transition) of the transition number matrix in the subject sequence. The “input pattern” is a pattern obtained by averaging and normalizing the sensor pattern (which will be described later). And the “cluster” is a specific number of patterns representative of all the input patterns, and are used for mapping an enormous types of input patterns into a specific number (30 in the present embodiment) of representative patterns.

[0079] The abnormality detection device 30 does not always require the distinction of whether the subject event belongs to the daily activity or the sample activity, and may process the stored and averaged data as the data in the daily activity. In other words, the data collecting unit 32 may store the collected data regardless of whether it is the data of the sample activity or the daily activity, and the Markov chain operating unit 33 may calculate above-mentioned characteristics by regarding the data stored for a fixed time period and averaged as the data in the daily activity and the individual (daily) data as the data in the sample activity.

[0080] The comparing unit 34 is a processing unit which is realized by a CPU or a memory that executes the data processing based on the control program. The comparing unit 34 has a function of comparing the daily activity and the sample activity, that is, the characteristics about the appearance frequency and the duration time of the state transition in Markov chain thereof, judging that unusual activity has occurred when the difference between them exceeds a predetermined threshold (or the similarity thereof is a predetermined threshold or less), and notifying the output unit 35 of the result, and has a distance calculating unit 34 a and a likelihood calculating unit 34 b.

[0081] The distance calculating unit 34 a calculates the distance in the appearance frequency of Markov chain based on the two transition number matrixes outputted from the Markov chain operating unit 33, that is, the transition number matrix of Markov chain obtained in the daily activity and the transition number matrix of Markov chain obtained in the sample activity. On the other hand, the likelihood calculating unit 34 b calculates the likelihood in the appearance frequency of Markov chain based on the transition number matrix of Markov chain obtained in the daily activity outputted from the Markov chain operating unit 33 and the cluster sequence obtained in the sample activity, or calculates the likelihood in the duration time distribution of Markov chain based on the duration time distribution of Markov chain obtained in the daily activity outputted from the Markov chain operating unit 33 and the cluster sequence obtained in the sample activity.

[0082] The output unit 35 is a CPU, a modem or the like that executes communication control based on the control program. Upon receipt of the notice from the distance calculating unit 34 a and the likelihood calculating unit 34 b in the comparing unit 34, that is, the notice of occurrence of unusual activity, the output unit 35 reports it to the pre-registered contact station (such as the monitoring center 50 and the cell phone 60) using the input unit 31 or the like via the communication network 40. For example, it makes a call using a cell phone, and sends a message indicating the occurrence of unusual activity to a pre-registered e-mail address.

[0083] The output unit 35 has option functions of, according to the prior setting by the input unit 31, (1) outputting not only the message indicating the occurrence of unusual activity but also the information indicating what type of unusual activity has occurred (such as the information indicating where the unusual activity has occurred, in appearance frequency and/or duration time distribution, or distance and/or likelihood, the value and the occurrence time thereof), (2) alarming or displaying the occurrence of unusual activity by an alarm bell or a display board mounted on the abnormality detection device 30, and (3) repeating the report until receiving the receipt data from the contact station.

[0084] Next, the operation of the abnormality detection system 10 configured as above will be explained in detail.

[0085]FIG. 4 shows values of sensor signals outputted from the small motion detection sensors 25 a˜25 c (a sequence of a sensor pattern) collected by the data collecting unit 32. More specifically, FIG. 4 shows sensor signals (“1” is indicated when a person exists and he is in small motion) outputted from the 20 small motion detection sensors 25 a˜25 c placed in the intelligent house 20, and the times when the sensor signals changed. Individual sensor pattern consists of 20-dimensional elements (“1” or “0”), and every time one element changes, a new sensor pattern is stored and accumulated in the sample activity data storage unit 32 a or the daily activity data storage unit 32 b.

[0086] These sensor patterns do not include the information on the arrangement of the sensors. “1” is indicated while each of the small motion detection sensors 25 a˜25 c is in the state of detecting small motion until detecting no motion, and “0” is indicated when it is in the state of detecting nothing. These binary information indicating whether the sensor (whose location is unknown) detects small motion or not and the information of the time when the event has occurred are the raw data inputted to the abnormality detection device 30.

[0087]FIG. 5 is a flowchart showing operations of the Markov chain operating unit 33. FIG. 5 shows an overall operation procedure of calculating the transition number matrix of Markov chain and the duration time distribution of Markov chain for the sequence of the sensor pattern shown in FIG. 4. Each step of this procedure will be explained in detail in the order of steps.

[0088] The Markov chain operating unit 33 first reads out the sequence of the sensor pattern from the sample activity data storage unit 32 a or the daily activity data storage unit 32 b depending upon the operation mode (monitoring mode and/or reference data creating mode), averages and normalizes it, and then transforms the sensor pattern into an input pattern.

[0089] More specifically, the sensor pattern b^(t) which is to be processed includes N steps (b^(t) is a sequence of b when t=1, . . . N). Since the information from the small motion detection sensors 25 a˜25 c sometimes occurs at intervals of several seconds, it is highly possible that appropriate weights cannot be assigned to human activities if such information is used as an input event as it is. Therefore, the obtained sensor pattern b^(t) is averaged in the time domain. Specifically, as shown in the following expression 1, the sensor pattern is divided into time windows, performed weighted addition respectively according to a Gaussian function, and then normalized.

[0090] sensor pattern: b^(t) (J-dimensional vector, t=1, . . . ,N)

[0091] averaged pattern: x^(t) (J-dimensional vector, t=1, . . . ,N)

[0092] variance: σ

[0093] data size: N

[0094] time window size: K $x^{\prime \quad t}:={\sum\limits_{k = 0}^{K - 1}{\left( {\frac{1}{\sqrt{2\pi}\sigma}^{- \frac{k^{2}}{2\sigma^{2}}}} \right) \cdot b^{t - k}}}$ $x^{t}:=\frac{x^{\prime \quad t}}{x^{\prime \quad t}}$

[0095] The detailed calculation procedure for averaging as shown in the first equation of the above expression is shown in the flowchart in FIG. 6, and the detailed calculation procedure for normalizing as shown in the second equation of the expression is shown in the flowchart in FIG. 7. In these figures, b^(t) and X^(t) are the data having 20-dimentionsional array elements and the duration time of the patterns thereof (b^(t·time, X) ^(t·time).)

[0096] Next, the Markov chain operating unit 33 vector-quantizes the obtained input pattern X^(t) so as to specify the representative input pattern (I pieces of clusters wi) as shown in the following expression.

[0097] cluster: w_(i) (J-dimensional vector, i=1, . . . ,I)

[0098] number of clusters: I

[0099] number of learnings: L[epoch] w_(c)^(t + 1) := w_(c)^(t) + ɛ(x^(t) − w_(c)^(t))

[0100] In the above expression, c:=arg max(wi·X^(t)) means that the vector that is most approximate to the vector X^(t) among the vectors wi(I=1, . . . ,I) (the vector that has the largest cosine of the angle between both vectors) is to be wc. The detailed calculation procedure for vector-quantization as shown in the expression is shown in the flowchart in FIG. 8. Here, the number of clusters I is 30 and the number of learnings L is 20,000.

[0101] As a result of the learning, the vector which is approximate to I pieces of clusters is calculated from the data space of the sensor patterns. The samples (clusters) as a result of it are shown as “x” marks in the plan view of the house in FIG. 9. Here, each cluster is represented as a center of gravity for the location of the small motion detection sensors 25 a˜25 c.

[0102] Next, as shown in the following expression, the Markov chain operating unit 33 obtains cluster sequence S^(t) corresponding to the input patterns X^(t) by inputting the input patterns X^(t) obtained in Step S10 into the clusters obtained in Step S11 so as to classify them. In other words, it converts the sequence of the input pattern into the sequence using specific pieces of clusters.

[0103] cluster sequence: s^(t) (s^(t)={1, . . . ,I},t=1, . . . ,N) $s^{t}:={\arg \quad {\max\limits_{i}\left( {w_{i} \cdot x^{t}} \right)}}$

[0104] The calculation procedure of the cluster sequence as shown in the above expression is shown in the flowchart in FIG. 10.

[0105] Next, the Markov chain operating unit 33 calculates D-dimensional Markov chain for the obtained cluster sequence S^(t). More specifically, it calculates the transition number matrix of Markov chain and the duration time distribution of Markov chain according to the following expression.

[0106] Calculation of transition number matrix

[0107] dimension of Markov chain: D

[0108] transition number matrix: M=└m_(i) _(⁰) _(, . . . ,i) _(^(D-1)) _(,i) _(^(D)) ┘

[0109] m_(i) _(⁰) _(, . . . ,i) _(^(D-1)) _(,i) _(^(D)) is the number of transitions from the state (i⁰ _(, . . . ,i) ^(D-1)) which appeared in s^(t) to the state i^(D).

[0110] The calculation procedure of the transition number matrix shown in the above expression is shown in the flowchart in FIG. 11. It is not regarded as transition when the clusters adjacent to each other in the cluster sequence are same, the processing for not counting such a case as a number of transitions is performed (Step S163).

[0111] Calculation of duration time distribution

[0112] Next, as for

m_(i) _(⁰) _(, . . . ,i) _(^(D-1)) _(,i) _(^(D))

[0113] the histogram obtained by dividing the time between the maximum value and the minimum value into L pieces as the duration time distribution indicating how long the state transition which appeared in the cluster sequences S^(t) continued

H_(i) _(⁰) _(, . . . ,i) _(^(D-1)) _(,i) _(^(D))

[0114] is generated.

[0115] number of divisions of histogram for duration time distribution: L

[0116] duration time distribution: H_(i) _(⁰) _(, . . . ,i) _(^(D-1)) _(,i) _(^(D)) =└h^(l) _(i) _(⁰) _(, . . . ,i) _(^(D-1)) _(,i) _(^(D)) ┘ (l=1, . . . ,L

[0117] h_(i) _(⁰) ^(l) _(, . . . ,i) _(^(D-1)) _(,i) _(^(D)) is the /th element of histogram for the duration time of transition from the state (i⁰, . . . ,i^(D-1)) which appeared in s^(t) to the state i^(D), where the time between the maximum value and the minimum value is divided into L pieces.

[0118] The calculation procedure of the duration time distribution is shown in the flowcharts shown in FIG. 12 and FIG. 13. Here, after calculating the time required for each state transition (Steps S170 ˜S178 in FIG. 12), the total number of events belonging to each time interval obtained by dividing the time between the maximum time and the minimum time into K pieces is calculated (Steps S180 ˜S187 in FIG. 13). In FIG. 13, Bar(k)·To and Bar(k)·From show a range of each time interval, and Bar(k)·height shows the total number of events belonging to the time interval.

[0119] The transition number matrix and the duration time distribution for the daily activities which are obtained by the above procedure are used for judging the occurrence of unusual activity in the comparing unit 34 as the daily activity template (reference data).

[0120]FIG. 14 is a data flow diagram showing a calculation procedure of distance in the appearance frequency of Markov chain performed by the distance calculating unit 34 a of the comparing unit 34.

[0121] First, the distance calculating unit 34 a calculates the probability matrix M′_(daily) by normalizing the transition number matrix M_(daily) calculated as the daily activity template (Step S20), and calculates the probability matrix M′_(sample) by normalizing the transition number matrix M_(sample) calculated from the sample activity sequence (Step S21). This probability distribution represents the appearance frequency of Markov chain.

[0122] Next, the distance calculating unit 34 a calculates the distance between these two probability distributions (Euclid distance in the present embodiment) according to the following equation (Step S22).

[0123] error=D(M′_(daily), M′_(sample))

[0124] The detailed calculation procedure for normalization in Steps S20 and S21 is shown in the flowchart shown in FIG. 15, and the calculation procedure of Euclid distance in Step S22 is shown in the flowchart shown in FIG. 16.

[0125] The larger the value of the distance “error” calculated as above is, the farther from the daily activity the person's activity is. Therefore, the comparing unit 34 judges whether this value “error” exceeds a predetermined threshold or not, and when it exceeds, notifies the output unit 35 of it.

[0126] In this comparison using the distance of the appearance frequency, the ratios of the numbers of activities are compared, and therefore, it is useful as rough observation of living activities. For example, if a person who has a meal after he gets home on a daily basis does not have it, the sequence where the activity of having meals is described drops out. This difference emerges noticeably in the transition number matrix, and the distances of the probability distributions that are the normalized transition number matrixes can be compared. However, since the time period dependence (the time range to be chosen) in calculating the transition number matrix is significant, the abnormality detection device 30 calculates the transition number matrix according to the time range specified by the operator who is designated via the input unit 31.

[0127]FIG. 17 is a data flow diagram showing a calculation procedure of likelihood in the appearance frequency of Markov chain performed by the likelihood calculating unit 34 b of the comparing unit 34.

[0128] First, the likelihood calculating unit 34 b calculates the probability distribution M′_(daily)

M′_(daily)=└m′_(i) _(⁰) _(, . . . ,i) _(^(D-1)) _(,i) _(^(D)) ┘

[0129] by normalizing the transition number matrix M_(daily) calculated as the daily activity template (Step S30).

[0130] Next, the likelihood calculating unit 34 b obtains the average log likelihood of the sample activity sequence (cluster sequence) s1, . . . ,sN to be compared with the above probability distribution as shown in the following expression to make it a characteristic for abnormality detection (Step S31). ${error} = {\frac{1}{N - D}{\sum\limits_{t = 1}^{N - D}{\log \left( m_{s^{t},\ldots \quad,\quad s^{t + D - 1},\quad s^{t + D}}^{\prime} \right)}}}$

[0131] The calculation procedure of likelihood in the above Step S31 is shown in the flowchart in FIG. 18. In this figure, the calculation procedure of the log likelihood “LineGraphs(t)·likelihood” for the cluster sequence to be compared therewith is shown. By calculating the average of the log likelihoods, the average logic likelihood “error” is calculated as shown in the expression 9.

[0132] The lower the value of the likelihood calculated as above is, the farther from the daily activity the activity indicated by the activity sequence to be compared is. Therefore, the comparing unit 34 judges whether or not this likelihood “error” is a predetermined threshold or less, and when it is the threshold or less, it notifies the output unit 35 of it.

[0133] This comparison using the likelihood of the appearance frequency comes into play for detecting abnormal condition where an unusual activity occurs over and over. For example, if a person who does not cook on a daily basis appears in the kitchen many times, it is detected as an unusual activity.

[0134]FIG. 19 is a data flow diagram showing a calculation procedure of likelihood for the duration time distribution of Markov chain performed by the likelihood calculating unit 34 b of the comparing unit 34.

[0135] First, for the duration time distribution of Markov chain of the daily activity template,

H_(i) _(⁰) _(, . . . ,i) _(^(D-1)) _(,i) _(^(D))

[0136] the likelihood calculating unit 34 b obtains the probability distribution obtained by normalizing the above duration time distribution

[0137] probability distribution H′h_(i) _(⁰) ^(l) _(, . . . ,i) _(^(D-1)) _(,i) _(^(D)) =└h′^(l) _(i) _(⁰) _(, . . . ,i) _(^(D-1)) _(,i) _(^(D)) ┘

[0138] (Step S40). The likelihood calculating unit 34 b obtains the average log likelihood of the activity sequence (cluster sequence) s1, . . . , sN to be compared with this probability distribution as shown in the following expression, and makes it a characteristic for detecting abnormality (Step S41). ${error} = {\frac{1}{N - D}{\sum\limits_{t = 1}^{N - D}{\log \left( h_{s^{t},\ldots \quad,\quad s^{t + D - 1},\quad s^{t + D}}^{\prime \quad l} \right)}}}$

[0139] The calculation procedure of likelihood in the above Step S31 is shown in the flowchart shown in FIG. 20. In this figure, “Barposition(Time)” in Step S225 is a function for checking what number the time “Time” is in the “Barso( )” array. This figure shows the calculation procedure of the log likelihood “LineGraphs(t)·likelihood” for the cluster sequence to be compared therewith. By calculating these log likelihoods, the average log likelihood “error” is calculated as shown in the expression 12.

[0140] The lower the value of the likelihood calculated as above is, the farther from the daily activity the activity indicated by the activity sequence to be compared is. Therefore, the comparing unit 34 judges whether or not this likelihood “error” is a predetermined threshold or less, and when it is the threshold or less, the comparing unit 34 notifies the output unit 35 of it.

[0141] In this comparison method using the likelihood for the duration time distribution, a characteristic indicating how long each activity takes on a daily basis is seen. For example, the activity of “passing by the stairs” continues only a few seconds on a daily basis and the duration time distributes around values of a few seconds. If there is no movement for a few minutes (the Markov chain does not continue), it is detected as an unusual activity.

[0142] Next, the calculation of the daily activity template and the details and results of five experiments for unusual activity detection using the template will be explained below.

(EXPERIMENT 1) MORNING SCENARIO

[0143] Assuming the activity sequence where a subject person who lives in the intelligent house 20 wakes up in the morning and goes out, the subject conducted the following activities and the sensory data thereof was stored. Samples of the transition number matrixes of Markov chains are shown in FIG. 21 and FIG. 22. FIG. 21 shows the transition number matrix in the daily activity template, and FIG. 22 shows the transition number matrix in the unusual activity. In these figures, full lines and broken lines which connect the clusters (x marks) indicate the activities which made a lot of transitions, that is, the full lines indicate the transitions from the clusters with smaller numbers to those with larger numbers and the broken lines indicate the transition from the clusters with larger numbers to those with smaller numbers.

[0144] Example of daily activities:

[0145] wake up→wash up→restroom→Ago out: 6 times

[0146] wake up→wash up→go out: 6 times

[0147] Example of unusual activities:

[0148] wake up→go out: 2 times

[0149] wake up→restroom→go out: 2 times

[0150] wake up→wash up→restroom→go out: 1 time (different subject)

[0151] In this experiment, the transition number matrix of Markov chain obtained by averaging the sensor data stored as daily examples was used as a daily activity template. FIG. 23 shows the comparison results made between the daily and unusual patterns and the daily activity template for the daily Markov chains using the distance comparison of the appearance frequency of Markov chain.

[0152] In FIG. 23, the abscissa indicates the numbers of the experiments and the ordinate indicates the Euclid distances. The experiment numbers 1 a 1˜1 b 6 are daily examples, and the experiment numbers 1 c 1˜1 a 7 k are unusual patterns. The experiment number 1 a 7 k is the unusual pattern of the different subject.

[0153] This result apparently shows that the distances of the unusual patterns from the daily activity template are longer and the distances of the daily examples from the daily activity template are shorter. Therefore, comparison between this distance and a predetermined threshold allows detection of unusual activity occurrence.

(EXPERIMENT 2) EVENING SCENARIO

[0154] Assuming that a person who lives in the intelligent house 20 comes back home, has dinner, relaxes in the sofa and then goes to bed, a subject person conducted the following activities and the sensor data thereof was stored.

[0155] Examples of daily activities:

[0156] come in→prepare dinner→dinner→wash dishes→relax in sofa →go to bed: 14 times

[0157] (Restroom action is added, relaxing timing may be shifted sometimes. For example: come in→restroom→prepare dinner→dinner→wash dishes→relax in sofa→go to bed)

[0158] Examples of unusual activities:

[0159] come in→go to bed

[0160] come in→restroom→go to bed

[0161] come in→wash up→go to bed

[0162] come in→prepare dinner→go to bed

[0163] come in→relax in sofa→restroom→go to bed

[0164] come in Japanese room→prepare dinner→dinner→wash dishes→go to bed

[0165] In this experiment, the transition number matrix of Markov chain obtained by averaging the sensor data stored as daily examples was used as a daily activity template. FIG. 24 shows the comparison results made between the daily and unusual patterns and the daily activity template for the daily Markov chains using the distance comparison of the appearance frequency of Markov chain.

[0166] In FIG. 24, the experiment numbers 2 a 1˜2 s 2 are daily examples, and the experiment numbers 2 k˜2 j 2 k are unusual patterns. The experiment number 2 j 2 k is the unusual pattern of the different subject.

[0167] This result apparently shows that the distances of the unusual patterns from the daily activity template are longer and the distances of the daily examples from the daily activity template are shorter. Therefore, comparison between this distance and a predetermined threshold allows detection of unusual activity occurrence.

(EXPERIMENT 3) SUDDEN FAINT

[0168] Assuming that a person who lives in the intelligent house 20 suddenly faints near the stairs or in the living room, the activities conducted when a subject person w alked around as usual and suddenly fainted near the stairs or in the living room and did not move were used as a sample. The sensor data collected in the above experiment 1/2 was used as a daily activity template.

[0169]FIG. 25 shows the comparison results between this sample and the daily Markov chains using the likelihood to the duration time distribution. In FIG. 25, the abscissa indicates the time intervals and the ordinate indicates the likelihood in these time intervals.

[0170] In this experiment, when the subject fainted in the living room (graph 400 in the figure) or on the stairs (402 in the figure), the values of the likelihood are low. Therefore, comparison between this likelihood and a predetermined threshold allows detection of unusual activity occurrence.

(EXPERIMENT 4) APARTMENT STYLE INTELLIGENT HOUSE

[0171] In this experiment, an apartment style intelligent house was used. In this intelligent house, 15 small motion sensors are placed, differently from the two-story experimental house, and a subject person actually lived there for 14 days. Assuming that the subject conducts a similar activity at the same time of the day for his life data, the transition number matrix obtained by averaging the activity sequences for 14 days divided into 1-hour time period (the full lines and broken lines in the plan view of the house in FIG. 26) was used as a daily activity template.

[0172]FIG. 27 shows the comparison results between the activity sequences of each day and this daily activity template using the likelihood to the appearance frequency of Markov chains. In FIG. 27, the abscissa indicates the times of the day and the ordinate indicates the dates. Also, the black squares in this figure show that there is no or very few activity, and the checked and white squares show that the tighter the check is, the higher the likelihood is.

[0173] In this figure, let us focus attention on 19˜20 o'clock of February 12 where the likelihood is extremely low. Examining the Markov chain in this time period in detail, the subject went to the washroom often. The subject went to a beauty salon in the afternoon of that day according to her report, it is supposed that that is why she went the washroom often.

(EXPERIMENT 5) DAILY ACTION RECORD

[0174] The daily action record is data based on self reports of 11 subject persons for one month regarding their 5 activities: sleeping, restroom use, having meals, going out and bathing (FIG. 28). FIG. 28 is a table showing the details of their reports, that is, the time periods of each activity, on a time series basis.

[0175] The entries of the subjects are acquired via the input unit 31 of the abnormality detection device 30 and operated by the Markov chain operating unit 33, and then the above data is realized. It is considered that this entry data shows their daily activity seen from more temporally macroscopic viewpoint than the data obtained by the small motion detection sensors 25 a˜25 c. Since this activity sequence is not the data obtained by the small motion detection sensors 25 a˜25 c, it is processed not based on vector-quantization, but only based on Markov chains.

[0176]FIG. 29 shows the result of the unusual activity detection for one subject performed using the likelihood comparison between the appearance frequency of Markov chains and the daily activities. The dates and times in this figure are same as those in FIG. 28. The average of the transition matrixes for one month was used as a daily activity template.

[0177] As shown in this figure, the result shows that the likelihood in the morning of July 22 is low. The data corresponding to this day (FIG. 28) shows a lot of sleeping time periods. Although the subject usually conducts regular activities, he was in bed through the morning on that day. The subject reported that he was sick with a cold on that day.

[0178] As described above, the abnormality detection device 30 describes the individual daily life and detects abnormality of his life via the information obtained by the sensors placed in the intelligent house. The transition number matrix of Markov chain and the duration time distribution of Markov chain are used for calculation of a daily activity template. It is shown that a number of abnormal cases can be detected by comparing the unusual activities and the daily activities using the likelihood and the distance.

[0179] Next, a modification of the present embodiment will be explained.

[0180] In the above embodiment, one daily activity template is used as the measure for judging unusual activities. In the modification below, a plurality of daily activity templates by time interval are created in advance, and the best daily activity template is selected for every time interval from among these templates for comparison with the sample activity so as to be a measure for judging unusual activities. In the following, a global daily activity template in the above embodiment is called a “global template”, and a daily activity template by time interval which will be explained below is called a “local template”.

[0181] The significance of introducing the local template is as follows. It cannot be said that the daily life of a person who lives in the intelligent house 20 which is described using only one global daily activity template represents his daily life accurately because a plurality of activities are described using only one template. Human activities can be seen as a collection of many repeatable activities such as walking down the stairs, bathing, restroom use and going from the living room to the bedroom. So, assuming that human life is made up of a plurality of daily activities and these activities are triggered according to the situation, the sensor sequence is transformed into discrete repeatable time intervals so as to generate a local daily activity template by time interval (local template).

[0182] Using a method of image processing for automatic sampling of time intervals, the time intervals with high correlation between templates are sampled by matching them using a time window. Thereby, the daily life representation with temporal hierarchy can be constructed.

[0183] The basic structure of the abnormality detection system according to the present modification is same as that of the abnormality detection system 10 according to the above embodiment, but different in that the Markov chain operating unit includes a processing unit for the local template.

[0184]FIG. 30 is a functional block diagram showing a detailed structure of the Markov chain operating unit 36 of the abnormality detection system according to the modification of the present invention. This Markov chain operating unit 36 includes a local template unit 37 that performs processing on a local template in addition to the functions of the Markov chain operating unit 33 of the above embodiment.

[0185] The local template unit 37 includes a local template generating unit 37 a that is activated in the reference data creating mode, a local template storage unit 37 b that stores the result generated by the local template generating unit 37 a and a local template predicting unit 37 c that is activated in the monitoring mode.

[0186] The local template generating unit 37 a generates a plurality of local templates by a predetermined time interval based on the daily activity data sent from the data collecting unit 32, calculates the characteristics thereof (the predicted probability and duration time of the local templates), and stores the local templates and the characteristics thereof in the local template storage unit 37 b. The local template predicting unit 37 c predicts to select the best local template for the inputted sample activity data from among a plurality of potential local templates for every time interval based on the predicted probability and duration time of the local templates stored in the local template storage unit 37 b. The Markov chain operating unit 36 outputs for every time interval the information (the transition number matrix and the duration time distribution) same as that in the above embodiment to the comparing unit 34 using the local template selected by the local template predicting unit 37 c as a measure. The comparing unit 34 compares them to judge unusual activities.

[0187] The detailed processing performed by the local template generating unit 37 a and the local template predicting unit 37 c will be respectively explained using the following flowcharts.

[0188]FIG. 31 is a flowchart showing the operation of the local template generating unit 37 a. The local template generating unit 37 a resolves the daily activity data into activity sequences (Step S300).

[0189] Here, the local template generating unit 37 a resolves into activity sequences by matching templates using a time window, which is an application of an image processing method.

[0190] More specifically, the sequence resolution is performed according to the following algorithm. The sensor sequence obtained in the intelligent house 20 are the time when the sensor value changes and the sensor pattern at that time. In this case, as shown in the small motion sensor sequences in FIGS. 32A˜D, the sensor pattern is transformed into the sensor pattern sequence per second (the maximum resolution of time). FIGS. 32A˜D show the sensor sequences detected by four small motion detection sensors, which represent a morning scenario.

[0191] It is assumed that very similar time intervals are sampled respectively from two small motion sensor sequences, b₁(t₁,j), b₂(t₂,j) and t₁=(1 . . . N₁), t2=(1 . . . N₂), j=(1, . . . ,J). Here, t₁ and t₂ are times per second, and j is a sensor channel of each of the small motion detection sensors 25 a˜25 c. This sensor sequence is binary, and can be processed as a binary image. Then, template matching is performed using a time window.

[0192] The size W of the time window is defined, and b₁(t₁ . . . t₁+W,j), and b₂(t₂ . . . t₂+W, j) are correlated with each other. The correlation is represented by the following expression used for the template matching. ${R\left( {t_{1},t_{2}} \right)} = \frac{\sum\limits_{t = 1}^{W}{\sum\limits_{j = 1}^{J}{{b_{1}\left( {{t_{1} + t},j} \right)}{b_{2}\left( {{t_{1} + t},j} \right)}}}}{R^{\prime}\left( {t_{1},t_{2}} \right)}$ $\begin{matrix} {{R^{\prime}\left( {t_{1},t_{2}} \right)} = \quad \left\{ {\sum\limits_{t = 1}^{W}{\sum\limits_{j = 1}^{J}{{b_{1}\left( {{t_{1} + t},j} \right)}{b_{1}\left( {{t_{1} + t},j} \right)}}}} \right\}^{1/2}} \\ {\quad \left\{ {\sum\limits_{t = 1}^{W}{\sum\limits_{j = 1}^{J}{{b_{2}\left( {{t_{2} + t},j} \right)}{b_{2}\left( {{t_{2} + t},j} \right)}}}} \right\}^{1/2}} \end{matrix} \times$

[0193] This is a normalized correlation, and can be used as a distance measure between images. This correlation is calculated in every interval between t₁=(1, . . . ,N₁−W) and t₂=(1, . . . ,N₂=W), and the obtained correlation R is seen as an image of N₁−W×N₂−W, as shown in the correlations in FIGS. 33A˜J. FIGS. 33A˜J are diagrams showing examples of the correlation R(t₁, t₂), and calculation results of correlation of all the times for 4 types of the sensor sequences as shown in FIGS. 32A˜D.

[0194] The image obtained by this correlation is binarized with a threshold α according to the following expression. ${R_{bin}\left( {t_{1},t_{2}} \right)} = \left\{ \begin{matrix} {{1\quad {R\left( {t_{1},t_{2}} \right)}} > \alpha} \\ {{0\quad {R\left( {t_{1},t_{2}} \right)}} \leq \alpha} \end{matrix} \right.$

[0195] This binary image R_(bin)(t₁, t₂) is labeled, and the labeled image is represented by L_(v), (v=1, . . . ,V). The starting time (the minimum value of each axis) t^(v) _(1,min) and t^(v) _(2,min) and the ending time (the maximum value) t^(v) _(1,max) and t^(v) _(2,max) are sampled for each axis t₁ and t₂ of L_(v),(v=1, . . . ,V), and (t^(v) _(1,min), t^(v) _(1,max)) and (t^(v) _(2,min), t^(v) _(2,max)) at that time are treated as the same sequences that are similar intervals of the sequence.

[0196] FIGS. 34A˜J shows examples of binarization R_(bin)(t₁, t₂) and labeling L_(v),(v=1, . . . ,V) of the correlation R(t₁, t₂) and sampling of each time period thereof. The sequence intervals with the same pattern in the interval indicated by arrows are the same sequence intervals.

[0197] Since this shows the correlation between two intervals only, the similarity in every combination of labeled images V×V needs to be considered. So, considering (t^(v) _(1,min), t^(v) _(1,max)) and (t^(v) _(2,min), t^(v) _(2,max)) as a two-dimensional vector, it is compared with the interval sampled from the other labeled image to regards the similar one as a similar interval. Euclid distance is used as this measure. A threshold μ is specified, and it is determined that they are the same intervals when the distance between them is the threshold or less. This is represented by the following expression. As for two labeled images v₁ and v₂, when

D(v ₁ , v ₂)={square root}{square root over ((t _(min) ^(v) ^(₁) −t _(min) ^(v) ^(₂) )+(t _(max) ^(v) ^(₁) −t _(max) ^(v) ^(₂) ))}<μ

[0198] is effected, it is considered that

(t_(min) ^(v) ^(₁) ,t_(max) ^(v) ^(₁) ) and (t_(min) ^(v) ^(₂) ,t_(max) ^(v) ^(₂) )

[0199] are similar intervals. These same similar intervals are collected as sets and respectively represented by u_(g),(g=1, . . . ,G). Note that u_(g) are disjoint sets to each other because the same similar intervals are treated as one set.

u ₁={(t _(1,min) ^(v) ^(₂) ,t _(1,max) ^(v) ^(₂) ), (t _(2,min) ^(v) ^(₂) ,t _(2,max) ^(v) ^(₂) ) ex., (t _(1,min) ^(v) ^(₃) ,t _(1,max) ^(v) ^(₃) ),(t _(2,min) ^(v) ^(₃) ,t _(2,max) ^(v) ^(₃) )}

[0200] FIGS. 35A˜H and FIGS. 36A˜D show examples of sampling of similar intervals for a complicated sequence. FIGS. 35A˜H show examples of sampling of similar sequence intervals in 4 cases of morning scenario. FIGS. 35A˜D respectively indicate the sensor sequences of the case 1 (wake up→wash up→restroom→go out), the case 2 (wake up→wash up→restroom→go out), the case 3 (wake up →restroom→go out) and the case 4 (wake up→go out), and FIGS. 35E ˜H respectively indicate the examples of sampling of similar sequence intervals corresponding to the cases 1˜4.

[0201] FIGS. 36A˜D show examples of sampling of similar sequence intervals in 2 cases of evening scenarios. FIGS. 36A and B respectively indicate the sensor sequences of the case 1 (come in →light ON→relax in sofa→prepare dinner→dinner→wash dishes→restroom→light OFF→go to bed) and the case 2 (come in→light ON →restroom→relax in sofa→prepare dinner→dinner→wash dishes→relax in sofa→light OFF→go to bed), and FIGS. 36C and D respectively indicate the examples of sampling of similar sequence intervals corresponding to the cases 1˜2.

[0202] Note that the parameters used for the above processing are empirically W=10, α=0.7 and μ=10.

[0203] Next, the local template generating unit 37 a averages the local templates in the similar sequence intervals obtained in the above Step S300 (Step S301 in FIG. 31). In other words, it averages local templates using the similar interval obtained in Step S300 as one local template.

[0204] More specifically, the transition number sequence from the state j to the state i in the cluster sequences S^(t) _(v)(vεu_(g)) obtained by clustering each similar interval is defined by the transition number matrix F^(v)=[f^(v) _(ij)].

[0205] Using this, the transition probability based on Markov chain representing a local template is calculated according to the following expression. $M^{g} = {\left\lbrack m_{ij}^{g} \right\rbrack = \left\lbrack \frac{\sum\limits_{{v \in u_{g}}\quad}f_{ij}^{v}}{\sum\limits_{i,j}{\sum\limits_{{v \in u_{g}}\quad}f_{ij}^{v}}} \right\rbrack}$

[0206] Then, the local template generating unit 37 a calculates the Markov chain between the local templates obtained in Step S301 (Step S302 in FIG. 31).

[0207] More specifically, when the local template is obtained in Step S301, the local template generating unit 37 a calculates the template transition probability M^(template) where the transition between the templates is regarded as a Markov chain. By transforming the template transition probability into a prior distribution to be a local template prediction distribution, the occurrence probability of the local template which is likely to occur at present can be predicted based on the appearance results of the local templates in the past sequences.

[0208] The appearance sequence of the local templates is z^(p) and the appearance time are y^(p). P indicates the times when the local templates changed from the current time, and is obtained by transforming the times t into discrete data according to the appearance order. The number of possible states is the number G of all the local templates, and the number of G×G combinations of the state transitions can be defined. Here, the number of transitions from the state j to the state i in z^(p) is defined by the transition number matrix F^(template)=[f_(ij) ^(template)]. When it is normalized into probability, the following expression is obtained. $M^{template} = {\left\lbrack m_{ij}^{template} \right\rbrack = \left\lbrack \frac{f_{ij}^{template}}{\sum\limits_{i,j}f_{ij}^{template}} \right\rbrack}$ ${\sum\limits_{i,j}m_{ij}^{template}} = 1$

[0209] This transition probability M^(template) is transformed into the predicted probability

M_(predict) ^(template)

[0210] by normalizing M^(template) when the state j at the time p−1 is fixed. This predicted probability is obtained by the same calculation as the prior probability in the state transition between the local templates. $M_{predict}^{template} = {\left\lbrack m_{{predict},\quad {ij}}^{template} \right\rbrack = \left\lbrack \frac{m_{ij}^{template}}{\sum\limits_{i}m_{ij}^{template}} \right\rbrack}$

[0211] Also, for every chain between the local templates, the average duration time thereof q_(ij) ^(template) is calculated. $q_{ij}^{template} = {\frac{1}{f_{ij}^{template}}{\sum\limits_{\{{{{p|y^{p + 1}} = i},{y^{p} = j}}\}}\quad \left( {y^{p + 2} - y^{p}} \right)}}$

[0212] Note that the ending time of the transition of p, p+1 is the appearance time of p+2.

[0213] Finally, the local template generating unit 37 a stores the predicted probabilities and the duration times of a plurality of the local templates obtained in Step S300 and the local templates obtained in Step S302 in the local template storage unit 37 b (Step S303 in FIG. 31).

[0214]FIG. 37 is a flowchart showing the operation of the local template predicting unit 37 c. When the local template storage unit 37 b actually uses the local template prediction model according to Markov chains for the sample activities, it identifies the appearance of the local template j using the distance between the transition probabilities of Markov chains, and predicts the next template based on the identified reliability.

[0215] More specifically, the local template predicting unit 37 c first makes the sequence b^(t) of the inputted sample activity into a cluster in the same manner as Steps S10˜S12 in the flowchart in FIG. 5 so as to obtain the cluster sequence s^(t) (Step S310).

[0216] Then, as for the cluster sequence s^(t) obtained by making the inputted sequence b^(t) into a cluster, the local template predicting unit 37 c calculates the Markov chain M(t, t−q_(ij)) in the time interval (t, t−q_(ij)) of the same length as each local template chain duration time from the current time (Step S311).

[0217] Next, the local template predicting unit 37 c calculates this Markov chain and the distance (cosine) of Markov chain Mg of each local template so as to be the reliability λ^(g)(t) of the local template g at the time t (Step S312).

λ^(g)(t)=vector(M′(t,t−q _(ij)))·vector(M′ ^(g))′ ${M^{\prime}\left( {t,{t - q_{ij}}} \right)} = \frac{M\left( {t,{t - q_{ij}}} \right)}{\left. ||{M\left( {t,{t - q_{ij}}} \right)} \right.||}$ $M^{\prime g} = \frac{M^{g}}{\left. ||M^{g} \right.||}$

[0218] Note that “vector” represents a function for transforming a matrix into a column vector.

[0219] Finally, the local template predicting unit 37 c predicts the current template φεG for the reliability obtained in the above step from the following predicted probability of Markov chain.

M_(predict) ^(template)

[0220] Based on the above, the current Markov chain M^(current)(t) is defined as follows. $\begin{matrix} {{M^{current}(t)} = {\left\lbrack {m_{ij}^{current}(t)} \right\rbrack = \left\lbrack \frac{\sum\limits_{g\varepsilon G}{{a^{g}(t)} \cdot m_{ij}^{g}}}{\sum\limits_{i,j}{\sum\limits_{g\varepsilon G}{{a^{g}(t)} \cdot m_{ij}^{g}}}} \right\rbrack}} \\ {{a^{g}(t)} = {\underset{p = 1}{\max\limits^{t}}{\max\limits_{\phi \quad {\varepsilon G}}\left( {{\lambda^{g}(p)} \cdot m_{{predict},g,\phi}^{template}} \right)}}} \end{matrix}$

[0221] Based on the above, the Markov chain operating unit 36 performs sequential calculation of M^(current)(t) when it uses the local template as a judgment measure.

[0222] FIGS. 38A˜H and FIGS. 39A˜H are diagrams showing examples of various local templates determined as above. FIGS. 38A˜H are diagrams showing examples of reliability λ^(g)(t) and predicted templates M^(current)(t). FIGS. 38A˜D respectively show the reliability λ^(g)(t) in the above-mentioned 4 cases, and FIGS. 38E˜H respectively show the local templates M^(current)(t) predicted based on the reliability shown in FIGS. 38A˜D.

[0223] FIGS. 39A˜H are diagrams showing examples of inputted cluster sequences s^(t) and the predicted templates M^(current)(t). FIGS. 39A˜D respectively show the cluster sequences s^(t) in the above-mentioned 4 cases, and FIGS. 39E˜H respectively show the local templates M^(current)(t) predicted for the cluster sequences s^(t) shown in FIGS. 38A˜D.

[0224] The Markov chain used for predicting local templates can be extended as follows by combining it with the prediction of cases.

[0225] The transition matrix of Markov chain between the templates for each sequence s^(t) _(v)(v=1, . . . ,V) which was used for generating M^(template) is calculated to be F_(v) ^(template). Then weighted average of the predicted value of each case is obtained for this Markov chain. $F^{template} = {\sum\limits_{v}{\xi^{t} \cdot F_{v}^{template}}}$

 ξ′=vector(M′(t))·vector(M′ ^(g))′ ${M^{\prime}(t)} = \frac{M(t)}{\left. ||{M(t)} \right.||}$ $M^{\prime g} = {\frac{M^{g}}{\left. ||M^{g} \right.||}.}$

[0226] Based on the above, the prediction distribution can be generated.

[0227] Also, there are the following two measures for deleting unnecessary local templates.

[0228] As the first measure, on the assumption that the local template obtained from each sensor sequence appears with regular frequency, the variance of number of the appearances in the sensor sequence of each local template is calculated. The local template with a larger variance value is deleted. Based on this noise deletion measure, the sequence whose sampled time interval is very short and is considered not to be appropriately related to the human-perceptible activity unit is deleted. However, since a certain level of dispersion of appearance frequency must be accepted, it is preferable to choose a modest threshold.

[0229] As the second measure, a measure

similar to AIC is preset to search for the required local template according to the SA method.

[0230]

=−(average likelihood of all sequences) +v·(number of templates)

[0231] (wherein v is a coefficient and empirically set to be 0.05.)

[0232] A combination of local templates that minimizes

is searched according to the SA method. Based on this noise deletion measure, the templates are chosen so that the entropy is as low as possible and the sequence is represented with as few templates as possible. FIGS. 40A˜F show examples of templates chosen. FIGS. 40A˜F show sequences for sampling 6 types of the local templates, FIGS. 40G˜I show local templates before selecting the local templates respectively corresponding to FIGS. 40A˜F, and FIGS. 40M ˜R show local templates after selecting the local templates respectively corresponding to FIGS. 40A˜F.

[0233] Next, the advantages of using the local template determined as above as an activity standard will be explained based on the comparison with the global template, using the following 3 experiments.

(EXPERIMENT 1) MORNING SCENARIO

[0234] This experiment is same as Experiment 1 (Morning scenario) in the above embodiment. So, the transition matrix of Markov chain in the daily example used for leaning and evaluation is same as that shown in FIG. 21, while the transition matrix of Markov chain in the unusual pattern is same as that shown in FIG. 22.

[0235] As for the cases when the global template is used for the daily activities and the unusual activities and the cases when the local template is used for both activities, when the average log likelihoods Err^(likelihood) of the activity sequences to the transition probability of the daily activity template are compared, the result shown in FIG. 41 was obtained.

[0236] In FIG. 41, the abscissa indicates experimental samples and the ordinate indicates the likelihood Err^(likelihood) of each sample to the daily activity template. This result shows that the likelihood becomes lower when using the local template in addition to the global template and the sensitivity to the unusual activity becomes higher (i.e., the likelihood for the unusual activity is lower than that for the daily activity). In other words, the result shows that the ability of detecting unusual activity is low when using only the global template but the unusual activity can be detected when using the local template.

(EXPERIMENT 2) EVENING SCENARIO

[0237] This experiment is same as Experiment 2 (Evening scenario) in the above embodiment. Assuming the activity sequence where a person who lives in the intelligent house 20 comes back home, has dinner, relaxes and then goes to bed, a subject person conducted activities (14 times of daily pattern and 1 time of 6 types of unusual patterns respectively) in Experiment 2 in the above embodiment and the sensory data thereof was stored.

[0238] As for the cases when the global template is used for the daily activities and the unusual activities and the cases when the local template is used for both activities, when the average log likelihoods Err^(likelihood) of the activity sequences to the transition probability of the daily activity template are compared, the result shown in FIG. 42 was obtained.

[0239] This result shows that the likelihood becomes lower when using the local template in addition to the global template and the sensitivity to the unusual activity becomes higher (i.e., the likelihood for the unusual activity is lower than that for the daily activity).

(EXPERIMENT 1/2) MORNING AND EVENING SCENARIOS

[0240] The experiment of the unusual activity detection (Err^(likelihood)) was carried out for the scenarios in the above Experiment 1 (Morning scenario) and Experiment 2 (Evening scenario). FIG. 43 shows the result thereof.

[0241] When the global daily activity templates are used, the above two cases are mixed and a method such as the distance Err^(distance) between the transition probabilities of Markov chains cannot be used. However, by predicting the local templates for situation dependency, an unusual activity can be detected stably even in the mixed cases.

[0242] As described above, in the present embodiment and the modification thereof, a personal daily life is described and unusual activities are detected based on the sensor information in the intelligent house 20. As characteristic amount for the daily activity template, the transition probability based on a Markov chain model and the duration time distribution of Markov chain events are used. Also, comparisons between the unusual activities and the daily activities using likelihood and distance measures show that some unusual patterns can be detected. As a result, it is suggested that the human living activities in a house like the intelligent house 20 can be represented by a model using probability distribution. Since image processing can be applied to sampling of local templates, it is also suggested that a lot of existing image processing technologies can be applied thereto.

[0243] The local template used for describing human activities can be seen as a chunk in cognitive psychology, and a model for triggering the chunk using prediction between chunks is used as an application. This can also be regarded as a cognitive human activity description model.

[0244] In the hierarchical structure of sequences based on the local templates used in the modification, the inclusion relation of activities is not clearly shown. As a method for solving it, statistical processing of simultaneous eventual probability between the local templates, starting time, ending time and others may be used, but it probably requires an enormous amount of experimental data.

[0245] There are two very significant problems in modeling probability in a housing environment: data specific to a user cannot always be obtained until services are provided, and frequency differences of activities are large. Modeling using a Bayesian approach is desired as a method for solving these problems. A stochastic graph representation like Bayesiannet (Breese, J. S., Construction of belief and decision networks, Computational Intelligence, Vol. 8, No. 4, pp. 624-647, 1992) seems to be effective because it is very applicable to introduction of prior knowledge and sequence resolution into an activity unit which is performed in the present embodiment. If a prior template based on a questionnaire about a user's occupation, age, life pattern is generated and used as a prior knowledge, the customizing period of the daily activity template which requires sampling of data specific to the user may be shortened. For construction of Bayesiannet, a method has been proposed such as a search (Suzuki J., Learning Bayesian Belief Networks Based on the Minimum Description Length Principle: Basic Properties. IEICE Trans. Fundamentals, Col. E82-A, No. 9, 1999 and G. F. Cooper and E. Herskovits, A Bayesian Method for Constructing Bayesian Belief Networks from Database, in Uncertainty in Artificial Intelligence' 91, UCLA, CA, pp. 86-94, 1991) using MDL standard (J. Rissanen, Stochastic Complexity and Modeling, The Annals of Statistics, Vol. 14, No. 3, pp. 1080-1100, 1986). Also, as a tool, BAYONET (Motomura and Hara, Stochastic System with Learning Mechanism from Database: BAYONET, The 12th Japanese Society for Artificial Intelligence, 1998) has been proposed. Bayesiannet, which must generally be a non-circular directed graph, may present a problem in representing a person's life. As a solution for this problem, a graph transformation using a clustering algorithm such as JunctionTree is under study. Also, in order to deal with a sequence, DBN (Dynamic Belief Network) for modeling a system having time delay (state transition) has been proposed (Dean, T. and Kanazawa, K., A Model for Reasoning about Persistence and Causation, Computational Intelligence, Vol. 5, No. 3, pp. 142-150, 1989). As mentioned above, the infrastructure for applying Bayesiannet to the intelligent house 20 has been established, so the future development thereof seems to be promising.

[0246] Since an applicable abnormality detection method depends upon abnormality to be detected under the present situation, an algorithm for automatically selecting an abnormality detection module if necessary needs to be developed. There still is a problem of setting a threshold for detection. It would be probably necessary to develop an application for high-level approach to a resident in the long term. Since it is believed that higher-level understanding of human beings would be necessary in view of this problem, cognitive knowledge obtained from the intelligent house 20 and interactivity between living things and study thereof are expected.

[0247] The abnormality detection system according to the present invention has been explained as above based on the embodiment and the modification thereof, but the present invention is not limited to these embodiment and the modification thereof.

[0248] For example, sensor patterns indicating human activities are not limited to the sensor signals outputted from the small motion detection sensors 25 a˜25 c. As shown in the abnormality detection system 110 in FIG. 44, sensor signals from a CO2 sensor 170 a, an open/close sensor 170 b, a small motion detection sensor 170 c, a passage detection sensor 170 d and so on and signals indicating operational states (power ON/OFF, status/positions of various switches, and so on) from various electrical household appliances such as an air conditioner 180 a, a refrigerator 180 b, a ventilating fan 180 c and a microwave oven 180 d may be used.

[0249] When a house has a home network 290 to which various sensors 270 and various electrical household appliances 280 are connected as shown in the abnormality detection system 210 in FIG. 45, it may be configured so that a controller 230 monitors the states of the sensors 270 and the electrical household appliances 280 via the home network 290 and sends the information indicating the detection results and operational states thereof to an abnormality announcement device or a home controller in the house via a communication network so as to announce and present the information, or the controller sends the information to the monitoring center and the abnormality detection device placed in the monitoring center detects the occurrence of abnormality in the house. In other words, the abnormality detection system of the present invention may be configured for decentralized processing so that the intelligent house collects data and the monitoring center analyzes and judges the occurrence of abnormality.

[0250] Also, as shown in the abnormality detection system 310 in FIG. 46, it may be configured so that an abnormality reporting device 335 is connected to a home network 390 separately and independently from the controller 330 of the home network 390. In other words, it may configured so that the abnormality reporting device 335 placed in the house collects the sensor signals from the small motion detection sensors 325 a˜325 c via the home network 390, performs the same processing as the abnormality detection device 30 of the above-mentioned embodiment, and reports the abnormality to the monitoring center 50 or the cell phone 60 via the communication network 40 when detecting it.

[0251] Furthermore, in the present embodiment, a fixed value such as data pre-stored in the daily activity data storage unit 32 b and the average value of one-month data is used as a reference value (daily activity template) for judging the occurrence of an unusual activity. However, a method may be used for updating (having it learn) the data based on the report from the resident of the intelligent house 20 (data with a teacher). For example, when a day ends, whether the day is usual or not is indicated in the input unit 31, and the Markov chain operating unit 33 updates the daily activity template regularly based on the report from the input unit 31. For example, when the transition number matrix and the duration time distribution in the latest 100 days which are reported as “daily” are to be operated, the Markov chain operating unit 33 calculates the average value of the days with more heavily weighting the later days and updates it as a new daily activity template. Then it judges the occurrence of unusual activity using the updated daily activity template as a reference value. Accordingly, the abnormality detection system which follows the secular changes in the human activity mode is realized.

[0252] As apparent from the above explanation, the abnormality detection device according to the present invention is a device that detects occurrence of abnormality in an event under observation, comprising: an input unit operable to acquire a sequence of an input pattern that is data which depends upon the event; a transition analyzing unit operable to analyze a characteristic of a transition in the acquired sequence of the input pattern; a comparing unit operable to compare the analyzed characteristic of the transition with a predetermined reference value, and judge that abnormality has occurred in the event when the characteristic and the reference value are not approximate to each other within a predetermined range; and an output unit operable to output occurrence of abnormality when the comparing unit judges that abnormality has occurred.

[0253] Accordingly, since the occurrence of abnormality is detected based on the sequence of the input pattern which depends upon the changes of the events to be monitored, it can be detected without depending upon the type and number of the event or the space (such as a place and a time period) for monitoring the event. More specifically, when monitoring the human activities, unusual activities in various manners can be detected without depending upon the time period or the place, differently from the detection of unusual activity only based on the time when the restroom is not used.

[0254] The input pattern which is to be the criterion is not determined fixedly and empirically, but the actually measured data in the normal state or the average of the data obtained for a certain period may be used. Thereby, abnormality can be accurately judged reflecting the individual characteristics (such as the individual differences in the normal state) of the person or the thing to be monitored.

[0255] Furthermore, by applying the Markov process to the sequence of the input pattern, the occurrence of abnormality may be judged based on the distance and the likelihood in the appearance frequency of Markov chain or the likelihood to the duration time distribution of Markov chain. Thereby, abnormality is monitored from multiple viewpoints such as the occurrence frequency (number of occurrence) and duration time of the event, differently from the abnormality judgment only based on the time when the restroom is not used, and various types of abnormality such as abnormal frequency and abnormal duration time can be detected with high accuracy.

[0256] As described above, the present invention enables to detect abnormality in the life of an elderly person who lives alone. Therefore, the practical value thereof is extremely high in the present society which is going rapidly toward aging and nuclear family. 

What is claimed is:
 1. An abnormality detection device that detects occurrence of abnormality in an event under observation, comprising: an input unit operable to acquire a sequence of an input pattern that is data which depends upon the event; a transition analyzing unit operable to analyze a characteristic of a transition in the acquired sequence of the input pattern; a comparing unit operable to compare the analyzed characteristic of the transition with a predetermined reference value, and judge that abnormality has occurred in the event when the characteristic and the reference value are not approximate to each other within a predetermined range; and an output unit operable to output occurrence of abnormality when the comparing unit judges that abnormality has occurred.
 2. The abnormality detection device according to claim 1, wherein the transition analyzing unit calculates appearance frequency for each type of the transition in the sequence of the input pattern as the analysis, and the comparing unit judges whether abnormality has occurred or not by comparing the input pattern acquired by the input unit with a reference input pattern as for the appearance frequency.
 3. The abnormality detection device according to claim 2, wherein the transition analyzing unit calculates probability distribution for each type of the transition as the appearance frequency, and the comparing unit calculates a distance between the input pattern acquired by the input unit and a reference input pattern as for the probability distribution, and judges that abnormality has occurred when the distance exceeds a predetermined value.
 4. The abnormality detection device according to claim 2, wherein the transition analyzing unit calculates probability distribution for each type of the transition as the appearance frequency, and the comparing unit calculates a likelihood of the input pattern acquired by the input unit to the reference input pattern as for the probability distribution, and judges that abnormality has occurred when the likelihood is a predetermined value or less.
 5. The abnormality detection device according to claim 1, wherein the transition analyzing unit calculates a duration time distribution that is a distribution of time required for the transition in the sequence of the input pattern as the analysis, and the comparing unit judges whether abnormality has occurred or not by comparing the input pattern acquired by the input unit with a reference input pattern as for the duration time distribution.
 6. The abnormality detection device according to claim 5, wherein the comparing unit calculates a likelihood of the input pattern acquired by the input unit to the reference input pattern as for the duration time distribution, and judges that abnormality has occurred when the likelihood is a predetermined value or less.
 7. The abnormality detection device according to claim 1 further comprising a clustering unit operable to transform the sequence of the input pattern acquired by the input unit into a cluster sequence that is a predetermined representative input pattern, wherein the transition analyzing unit analyzes the characteristic in the cluster sequence which is transformed by the clustering unit.
 8. The abnormality detection device according to claim 7, wherein the clustering unit obtains the cluster sequence after averaging and normalizing the sequence of the input pattern acquired by the input unit.
 9. The abnormality detection device according to claim 8, wherein the clustering unit obtains the cluster sequence after specifying all the clusters by vector-quantizing the sequence of the input pattern obtained by the normalization.
 10. The abnormality detection device according to claim 1, wherein the reference value is a value, as for a reference event, obtained by acquiring the sequence of the input pattern in advance through the input unit and analyzing the characteristic through the transition analyzing unit.
 11. The abnormality detection device according to claim 10, wherein the reference value is a value determined based on data given by an operator which is learned with a teacher and the input pattern which is acquired repeatedly by the input unit.
 12. The abnormality detection device according to claim 1, wherein the reference value is a value specified by an operator.
 13. The abnormality detection device according to claim 1, wherein the input unit includes a plurality of small motion detection sensors placed in a plurality of places in a house, and acquires, as the input pattern, a sensor pattern that is a combination of data indicating whether there is a small motion or not outputted from a plurality of the small motion detection sensors.
 14. The abnormality detection device according to claim 13, wherein the input unit assumes that a new sensor pattern has occurred every time the sensor pattern changes and acquires a sequence of an input pattern corresponding to a sequence of the sensor pattern.
 15. The abnormality detection device according to claim 14, wherein the input unit acquires the input pattern by averaging and normalizing the sensor pattern in the time domain.
 16. The abnormality detection device according to claim 1, wherein the input unit acquires the sequence of the input pattern from an operator.
 17. The abnormality detection device according to claim 1, wherein the output unit reports to a predetermined destination via a transmission channel that abnormality has occurred.
 18. The abnormality detection device according to claim 1, wherein the transition analyzing unit calculates a Markov chain in the sequence of the input pattern and analyzes a characteristic of the calculated Markov chain as the analysis.
 19. The abnormality detection device according to claim 18, wherein the transition analyzing unit calculates an appearance frequency of the Markov chain as the analysis, and the comparing unit judges whether abnormality has occurred or not by comparing the input pattern acquired by the input unit with a reference input pattern as for the appearance frequency.
 20. The abnormality detection device according to claim 19, wherein the comparing unit calculates an Euclid distance between the input pattern acquired by the input unit and the reference input pattern as for the appearance frequency of the Markov chain, and judges that abnormality has occurred when the Euclid distance exceeds a predetermined value.
 21. The abnormality detection device according to claim 20, wherein the comparing unit calculates a likelihood of the input pattern acquired by the input unit to the reference input pattern as for the appearance frequency of the Markov chain, and judges that abnormality has occurred when the likelihood is a predetermined value or less.
 22. The abnormality detection device according to claim 18, wherein the transition analyzing unit calculates a duration time distribution of the Markov chain as the analysis, and the comparing unit judges whether abnormality has occurred or not by comparing the input pattern acquired by the input unit with the reference input pattern as for the duration time distribution.
 23. The abnormality detection device according to claim 22, wherein the comparing unit calculates a likelihood of the input pattern acquired by the input unit to the reference input pattern as for the duration time distribution of the Markov chain, and judges that abnormality has occurred when the likelihood is a predetermined value or less.
 24. The abnormality detection device according to claim 1 wherein the reference value is a collection of local reference values in each time interval which satisfies a predetermined condition, and the comparing unit judges whether abnormality has occurred or not by comparing the characteristic with the local reference value in said each time interval.
 25. The abnormality detection device according to claim 24 further comprising: a local reference value generating unit operable to generate the local reference value in said each time interval by acquiring a plurality of the sequences of the input patterns through the input unit, resolving the sequences in each time interval where the acquired input patterns are similar to each other within a predetermined range and collecting the resolved sequences in said each time interval; and a local reference value selecting unit operable to predict and select an optimum reference value from among all the local reference values belonging to said each time interval based on the input patterns for abnormality detection, wherein the comparing unit judges whether abnormality has occurred or not by comparing the characteristic with the local reference value selected by the local reference value selecting unit in said each time interval.
 26. The abnormality detection device according to claim 25, wherein the local reference value generating unit samples a time interval with a high correlation according to template matching using a time window and resolves the sequence in said each sampled time interval.
 27. The abnormality detection device according to claim 25, wherein the local reference value selecting unit calculates probability distributions for every transition type in the sequences of the input patterns for abnormality detection and the input patterns of the reference event respectively, and predicts and selects the optimum local reference value based on a distance between the probability distributions.
 28. An abnormality detection system comprising: an abnormality detection device that is placed in a place where an event for abnormality detection occurs; a communication device that monitors occurrence of abnormality; and a transmission channel that connects the abnormality detection device and the communication device, wherein the abnormality detection device includes: an input unit operable to acquire a sequence of an input pattern that is data which depends upon the event; a transition analyzing unit operable to analyze a characteristic of a transition in the acquired sequence of the input pattern; a comparing unit operable to compare the analyzed characteristic of the transition with a predetermined reference value, and judge that abnormality has occurred in the event when the characteristic and the reference value are not approximate to each other within a predetermined range; and an output unit operable to output occurrence of abnormality when the comparing unit judges that abnormality has occurred, and the communication device includes: a receiving unit operable to receive a report from the abnormality detection device; and a showing unit operable to show an operator the receipt of the report when the receiving unit receives the report.
 29. An abnormality detection system comprising: an abnormality detection device that is placed in a place where an event for abnormality detection occurs; a communication device that monitors occurrence of abnormality; and a transmission channel that connects the abnormality detection, device and the communication device, wherein the abnormality detection device includes: an input unit operable to acquire a sequence of an input pattern that is data which depends upon the event; and a sending unit operable to send the acquired input pattern to the communication device via the transmission channel, and the communication device includes: a receiving unit operable to receive the input pattern sent from the abnormality detection device; a transition analyzing unit operable to analyze a characteristic of a transition in the acquired sequence of the input pattern; a comparing unit operable to compare the analyzed characteristic of the transition with a predetermined reference value, and judge that abnormality has occurred in the event when the characteristic and the reference value are not approximate to each other within a predetermined range; and a showing unit operable to show an operator the occurrence of abnormality when the comparing unit judges that abnormality has occurred.
 30. The abnormality detection system according to claim 28 or claim 29, wherein the input unit includes a plurality of small motion detection sensors placed in a plurality of places in a house, and acquires, as the input pattern, a sensor pattern that is a combination of data indicating whether there is a small motion or not outputted from a plurality of the small motion detection sensors.
 31. The abnormality detection system according to claim 28 or claim 29, wherein the input unit includes an operation state detection sensor that detects an operation state of equipment which is placed in the place, and acquires a sensor signal from the operation state detection sensor as the input pattern.
 32. The abnormality detection system according to claim 28 or claim 29, wherein the place is a house, the abnormality detection system further comprises: a home network that connects a plurality of electrical household appliances which are placed in the house; and a controller that controls a plurality of the electrical household appliances via the home network, and the input unit detects operation states of a plurality of the electrical household appliances via the home network and acquires the detected operation states as the input patterns.
 33. The abnormality detection system according to claim 28 or claim 29, wherein the place is a house, the abnormality detection system further comprises a home network that connects a plurality of electrical household appliances which are placed in the house, the abnormality detection device further includes a controller that controls a plurality of the electrical household appliances via the home network, and the input unit detects operation states of a plurality of the electrical household appliances via the home network and acquires the detected operation states as the input patterns.
 34. An abnormality detection method for detecting occurrence of abnormality in an event under observation, comprising: an inputting step for acquiring a sequence of an input pattern that is data which depends upon the event; a transition analyzing step for analyzing a characteristic of a transition in the acquired sequence of the input pattern; a comparing step for comparing the analyzed characteristic of the transition with a predetermined reference value, and judging that abnormality has occurred in the event when the characteristic and the reference value are not approximate to each other within a predetermined range; and an outputting step for outputting occurrence of abnormality when the comparing unit judges that abnormality has occurred.
 35. The abnormality detection method according to claim 34, wherein in the transition analyzing step, appearance frequency for each type of the transition in the sequence of the input pattern is calculated as the analysis, and in the comparing step, whether abnormality has occurred or not is judged by comparing the input pattern acquired in the inputting step with a reference input pattern as for the appearance frequency.
 36. The abnormality detection method according to claim 34, wherein in the transition analyzing step, a duration time distribution that is a distribution of time required for the transition in the sequence of the input pattern is calculated as the analysis, and in the comparing step, whether abnormality has occurred or not is judged by comparing the input pattern acquired in the inputting step with a reference input pattern as for the duration time distribution.
 37. The abnormality detection method according to claim 34 further comprising a clustering step for transforming the sequence of the input pattern acquired in the inputting step into a cluster sequence that is a predetermined representative input pattern, wherein in the transition analyzing step, the characteristic in the cluster sequence which is transformed in the clustering step is analyzed.
 38. The abnormality detection method according to claim 34, wherein the reference value is a value, for a reference event, obtained by acquiring the sequence of the input pattern in advance in the inputting step and analyzing the characteristic in the transition analyzing step.
 39. The abnormality detection method according to claim 34, wherein the reference value is a value specified by an operator.
 40. The abnormality detection method according to claim 34, wherein in the inputting step, a plurality of small motion detection sensors placed in a plurality of places in a house are included, and a sensor pattern that is a combination of data indicating whether there is a small motion or not outputted from a plurality of the small motion detection sensors is acquired as the input pattern.
 41. The abnormality detection method according to claim 34, wherein in the inputting step, the sequence of the input pattern is acquired from an operator.
 42. The abnormality detection method according to claim 34, wherein in the outputting step, it is reported to a predetermined destination via a transmission channel that abnormality has occurred.
 43. The abnormality detection method according to claim 34, wherein in the transition analyzing step, a Markov chain in the sequence of the input pattern is calculated and a characteristic of the calculated Markov chain is analyzed as the analysis.
 44. The abnormality detection method according to claim 34, wherein the reference value is a collection of local reference values in each time interval which satisfies a predetermined condition, and in the comparing step, whether abnormality has occurred or not is judged by comparing the characteristic with the local reference value in said each time interval.
 45. A program for detecting that abnormality has occurred in an event under observation, the program causing a computer to execute: an inputting step for acquiring a sequence of an input pattern that is data which depends upon the event; a transition analyzing step for analyzing a characteristic of a transition in the acquired sequence of the input pattern; a comparing step for comparing the analyzed characteristic of the transition with a predetermined reference value, and judging that abnormality has occurred in the event when the characteristic and the reference value are not approximate to each other within a predetermined range; and an outputting step for outputting occurrence of abnormality when the comparing unit judges that abnormality has occurred. 